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Abstract 

A system put in contact with a large heat bath normally thermalizes. This 
means that the state of the system p s {t) approaches an equilibrium state pf q , the 
latter depending only on macroscopic characteristics of the bath (e.g. temperature), 
but not on the initial state of the system. The above statement is the cornerstone 
of the equilibrium statistical mechanics; its validity and its domain of applicability 
are central questions in the studies of the foundations of statistical mechanics. In 
the present paper we concentrate on one aspect of thermalization, namely, on the 
system initial state independence (ISI) of pf q . A necessary condition for the system 
ISI is derived in the quantum framework. We use the derived condition to prove 
the absence of the system ISI in a specific class of models. Namely, we consider 
a single spin coupled to a large bath, the interaction term commuting with the 
bath self-Hamiltonian (but not with the system self-Hamiltonian) . Although the 
model under consideration is nontrivial enough to exhibit the decoherence and the 
approach to equilibrium, the derived necessary condition is not fulfilled and thus 

depends on the initial state of the spin. 

1 Introduction 

In the last decade a considerable progress in deriving the fundamentals of statistical 
physics from the first principles of quantum theory was achieved PQ-[H] (see [H] for a 
profound list of references to the related works). Most of the results which contribute 
to this progress were obtained in the following framework. Consider a quantum system 
(described by a Hilbert space S), which interacts with a quantum bath (described by a 
Hilbert space B). The bath is considered to be "much larger" than the system. In par- 
ticular, if S and B are finite-dimensional with the dimensions d$ and dg correspondingly, 
which is assumed in what follows, then ds "C dg. The composite system with Hilbert 
space H = S (g) B is considered to be closed and to evolve according to the Shroedinger 
equation with a Hamiltonian 

H = H S + H B + H SB , (1) 
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where H s and H B are self-Hamiltonians of the system and the bath correspondingly, and 
is an interaction Hamiltonian. Here and in what follows the usage of superscripts 
and subscripts S, B is believed to be self-explanatory. The state of the combined system 
W is described by a state vector f G H. The latter evolves as ^(i) = exp(— \Ht)^(fS). 
The states of the system S and the environment B are described by the reduced density 
matrices 

pP=tT B \*)(#\, p B = ti s \^){^\ (2) 

correspondingly 

The initial state of the composite system H is taken to be a product state: 

(0) = ^ e 5, $ G B. (3) 

Here and in what follows we use Greek letters ip and $ to denote the normalized 
vectors of unitary spaces "H, S and £> correspondingly. The product form of the initial 
state is natural when considering the approach to the thermal equilibrium. Usually all 
results obtained with the use of the product initial state assumption may be generalized 
to the case of an arbitrary pure initial state. Moreover, usually the results obtained for a 
pure state \1/ are straightforwardly generalized to the case when the state of the composite 
system TL is mixed and described by the density matrix p. 

What can be said about the long-time behavior of p s {t) in the case of a generic 
interaction H SB 1 This is a central questions of the equilibrium statistical mechanics. An 
intuitive answer is that p s approaches an equilibrium density matrix of some special 
(e.g. canonical) form. As was argued in [5J, on closer examination one expects that the 
system exhibits four distinct properties, which we refer to as thermalization properties. 
We formulate them below exploiting the product form of the initial state ([3]): 

1. Equilibration. By definition, a system equilibrates if p s (t) approaches a time- averaged 
density matrix p s and stays close to it most of the time. Defined in this way equi- 
libration does not imply neither any special form of p s , nor the independence of p s 
from initial conditions. 

2. Bath initial state independence (Bath ISI). This means that p s (almost) does not 
depend on the exact microstate of the bath, $. Rather p s should depend on some 
macroscopic characteristics of the state of the bath, which should be represented 
by functionals on B. The prime example of such characteristic is the bath inverse 
temperature (3 = /?[$]• 

3. System initial state independence (System ISI). This means that p s (almost) does 
not depend on ip. 

4. Boltzmann-Gibbs form of the equilibrium state: p s = Z~ x exp(— f3H s ). This prop- 
erty may be expected if the interaction H is in some sense "weak" compared to the 
system self-Hamiltonian H s (although "sufficiently strong" to ensure equilibration) 
and the initial state of the bath has a small energy uncertainty. 

The last three properties make sense only if the firsts one holds. The last property makes 
sense if also the properties (2) and (3) hold. Note the lack of the symmetry between the 
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definitions of the bath ISI and the system ISI. This asymmetry arises because the bath is 
assumed to be much larger than the system. 

The ultimate goal is to derive all four properties from the first principles of quantum 
theory under reasonable conditions. 

The present paper addresses the third property. To start with, we briefly review the 
main results concerning all four properties. The first key fact was discovered in PQ-[I]. 
It is based on the concentration of measure phenomenon, which is a striking feature of 
geometry in spaces of very high dimensions. Consider a linear subspace Hr C % with 
dimensionality ^> 1. Then for almost all states \I/ G Hr the reduced density matrix 
trg|\I>) is close to the averaged over Hr matrix: 

tr B |tf)(*|~(tr B |tf)<*| ) Hr . (4) 

Usually the trajectory ^(t) entirely lies in some "Hr, for example in the energy shell, which 
is spanned by the eigenvectors of H with the eigenvalues in some range [E, E + AE] . In 
this case it is natural to assume that p s (t) = tr B |#(t))(tf(*)| for almost every t is close to 
the "H^-average ( tr#|\l/)(\l/| )% R . Such assumption implies that the property (1) generally 
holds; the properties (2) and (3) also hold provided that all the considered initial states 
belong to Hr. Also, under certain assumptions, one may perform the averaging over %r 
explicitly to get the Boltzmann-Gibbs form of the averaged reduced density matrix, which 
addresses the property (4). 

Although the above arguments provide an important insight into the problem, they 
do not constitute the proofs of properties (1) - (4). The reason is that the trajectory 
ty(t) actually never completely fills any high- dimensional linear subspace 1-Lr. Rather it 
generically fills some torus [10J. The dimensionality of the latter depends on the rate of 
the entanglement between the bath and the system introduced by the interaction H SB . 
In fact, the entanglement appears to be of a primary importance in the problem involved. 
In case of non-interacting system and environment (H SB = 0), when the entanglement is 
completely absent, none of the properties (1) - (4) hold. They do not also hold in case of 
non-zero but very weak interaction, when the operator norm of H is much smaller than 
the typical energy level spacing of the environment. In the latter case the perturbation 
theory may be used to calculate the eigenstates and eigenvalues of H and to demonstrate 
that p (t) does not equilibrate. 

Substantial success in derivation of the properties (1), (2) and, partly, (3) was achieved 
in [6j. A general quantum-mechanical closed system divided in two parts was considered; 
the only requirement for the total Hamiltonian H was the nondegeneracy of energy gaps. 
First, it was proven that the equilibration property holds provided the initial state ^(0) 
is a superposition of a large number of eigenvectors of the total Hamiltonian H. Second, 
it was shown that the equilibrium density matrix p s is almost the same for almost all 
initial states <£> of environment which belong to a high-dimensional subspace Br C £>; 
those states of the bath which provide the exceptions from this rule form a subset in 
Br of an exponentially small measure. Thus, the bath ISI property was proven. Third, 
an inequality was derived which in certain circumstances (in particular, when d$ ^> 1) 
proves the system ISI. However, as was emphasized in [BJ, the problem of the system 
state independence is rather complicated. In particular, it was pointed out in [9] that 



3 



the inequality derived in [6] is not very restrictive when the dimensionality of the Hilbert 
space of the system S is low. 

Among other important advances one should mention the strong upper bounds on the 
speed of fluctuations around the equilibrium state derived in |5j [8] . 

In the present paper the system initial state independence problem is addressed, es- 
pecially in the case when ds is small. In particular, a single spin 1/2 is considered as a 
system (d$ = 2). The paper is organized as follows. In section [2] we introduce definitions 
and notations. In section [3] we quote and discuss a key theorem from [B] which allows 
to prove the bath ISI and, in some cases, the system ISI. In section H] we present our 
main results concerning the system ISI property. In section [5] a class of exactly solvable 
models is considered for which the derived necessary condition is not satisfied and thus 
the system ISI property does not hold. The outline of our results is presented in section 
El 



2 Definitions 

Let us start from introducing the setup, definitions and notations. 

Any finite system returns to the arbitrarily small vicinity of its initial state infinitely 
many times (although the recurrence time is normally very large). For this reason the 
limit lim p s (t) does not exist. Instead a time-averaged density matrix 

t— >OD 

i r* 

f = lim - / p s {t')dt'. (5) 



t—>oo t 







should be considered. Note that throughout the paper we use an overline to denote the 
time- averaging, and angle brackets to denote the averaging over normalized vectors from 
some subset of a Hilbert space with a uniform measure, the latter being defined in the 
end of the present section. 

To quantify the difference between two states p\ and p<i we use the distance \\P\ — fall, 
where 

\\p\ \ = tTy/f?. (6) 

The maximal value of this distance is 2. This is a physically meaningful definition, as it 
is discussed e.g. in |J1 US]- In particular, it is equal to the doubled maximum difference 
in probability for any outcome of any measurement performed on the two states |6J. 
The total Hamiltonian is assumed to be nondegenerate, 

d 

H = Y,E n \^n)(^n\, (7) 

n=l 

where \l/ n are the eigenvectors of H, and E n < E m for n < m. The initial sate of the 
composite system is a superposition of the eigenvectors with coefficients c n : 

d 

$(0) = ^^, (8) 

n=l 



1 Note that the above defined distance is the same as in [3] but is two times larger than in [B\ The 
accepted definition is natural when considering the distance between the states of a single spin, see below. 
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The time evolution of p s {t) reads 

^w = EE c ^ e_i(E "" £m) VL, (9) 

n=l m=l 

where the matrices 

Pnm = tr B |tt ft )<* m | (10) 

are introduced (not to be confused with matrix elements!). Evidently, these matrices 
encode the dynamics of the open system S, while the coefficients c n describe the initial 
conditions. The nondegeneracy of the Hamiltonian allows to obtain 

n=l 

where a short-form notation p% = p^ n is used. The time-averaged state of the system p s 
depends, in general, on the initial states of the system and the bath, ip and $ correspond- 
ingly, through the coefficients c n = (\l/ ri |V'^ > ) : P s = p s [ip&]- 

In the case when the system S is represented by a single spin, any p s may be param- 
eterized by a polarization vector p, 

p 5 = (l + p0 -)/2, p = tr 5 (p 5 <r), 0<|p|<l (12) 

The polarization vector belongs to a unit sphere which is known as the Bloch sphere. The 
length of a polarization vector equals 1 for a pure state and is less than 1 for a mixed 
state. The distance between two states pf and pf is simply the Euclidian distance in the 
Bloch sphere: 

Ilpf-Pfll = IP1-P2I <2. 

The scalar product of vectors p and p' is denoted as (p, p')- We define the following 
important polarization vectors: the initial polarization vector po = (ipl&l'ip) , polarization 

vectors which correspond to eigenstates of composite system p n = (^ n |<r|^ n ) and the 

d 

time-averaged polarization vector p = X) l c n| 2 Pn- 

71=1 

In order to introduce averages and to formulate propositions about states which are 
typical for some subspace %r e H, we need to define a uniform measure on Hr. Strictly 
speaking, pure states of a physical system are in one-to-one correspondence with one- 
dimensional linear subspaces of a Hilbert space or, equivalently, with rank one projectors 
|\E')(\I/|. Therefore actually one should consider the projective space %rP instead of the 
Hilbert space 1-Lr. It is possible to define a uniform measure on a projective space of pure 
quantum states through the Haar measure on a SU(cIr) group, taking into account that 
any pure state |^ r )(^ r | may be obtained from some fixed state |\E'o)(^ / o| by the unitary 
transformation (see e.g. j9] for the details). However, following [3l |6] we use a different, 
more explicit construction, which leads to the same result. Namely, let us choose an 
arbitrary basis {^i} in Hr and establish a map 

$f>i6 R 2dR : x 2 i-i = Re(^|^), x 2 i = Im(ty|tf). (13) 



5 



All normalized vectors from TLr are therefore in one-to-one correspondence with points of 
the 2dji — 1-dimensional unit sphere embedded in the 2 cLr- dimensional Euclidian space. 
Note however that a physical state corresponds not to a point but to a one- 

dimensional curve on the sphere, because of the overall phase ambiguity of Now to 
pick up a quantum state from %r (more precisely, from HrP) at random according 
to the uniform measure, we first pick up a vector x from a unit sphere according to 
the uniform measure on a sphere, and then construct the corresponding state |^ r )(^ r |. A 
thereby constructed measure does not depend on the choice of the basis {^i}- 

In the above paragraph we reminded a well-known fact that actually a pure physical 
state should be characterized by a projector (or by a vector \I/ "up to a phase 

factor"). Bearing this in mind, in what follows we use a common language and speak 
about "state vectors" \l/ and "state spaces" "H,"H_r, ... without further stipulations. 



3 Sufficient condition for the system initial state in- 
dependence 

The following theorem concerning the initial state independence was proven in [6] : 
Theorem 0. Consider the hamiltonian H with nondegenerate energy gaps, which means 
that Ek — Ei = E m — E n implies either k = I, m = n, or k = m, I = n. 
(i) Almost all initial states chosen from a large restricted subspace 1-Lr C H with the 
dimensionality d,R yield the same equilibrium state. In particular, 



(p S ) 



Hi 




with 



5 ^(^nKdRj-'Un^n) tr s (p s n ) 



< 1, 



(14) 



(15) 



n=l 



where Hn R is the projector onto 1-Lr% 

(ii) There are exponentially few states in Hr which yield a substantial distance between 
p s and (p s )-R R . In particular, for a random state \l/ G %r 



Pr 



(P S ) 



H, 



> 



ds5 



d 



R 



< 2e~ cdRe \ 



(16) 



where c = l/(187r 3 ). 

First we review how this theorem may be used to prove the bath ISI property [6]. Let 
us choose any system state ip and consider Hr as a tensor product of if) and some large 
Br C B with the dimensionality dR^> d$ Hr = if> <8> Br. Then one gets that for a fixed 
if) and vast majority of $ G Br the equilibrium state p s is close to the average (p s )b r - in 
other words, it is proven that for any fixed initial state if) of the system the equilibrium 
state p s depends on the initial state of the bath $ G Br extremely weakly. 



2 In general, the initial condition is not required to be of a product form in this theorem. 
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Note that the smallness of 5 is not required in the above proof; in fact one may safely 
take 6 = 1 and exploit the weaker bound in (j!4p . The bound with S was introduced in 
[5] in order to treat the system ISI problem. The latter appears to be more complicated 
compared to the previous one. Indeed, let us try to proceed analogously to what was done 
in the preceding paragraph. We fix some state of the bath $ and construct T-Lr = S £g> $ 
(according to the formulation of the system ISI property we should take the whole S 
instead of some subspace in S). Now, however, cIr = d s , and the weaker bound in (TTJJ 
appears to be useless. The stronger bound is useful provided s/5 is small. For this reason 
one may look at Theorem with Hr = S <g> $ as on the 

Sufficient conditions for the system ISI: If 

V / 5<1, (17) 

then the system ISI property holds. 

In [5] a case was considered when, firstly, the dimensionality of the system is large, 
yfds 3> 1, and, secondly, the eigenstates \l/ n are highly entangled (in particular, far from 
product), which implies that the purities of the density matrices are close to their 
minimal values: 

tr s (p s n ) 2 ~ l/ds. (18) 

In this case S ~ l/ds, the above sufficient condition is satisfied and the system ISI property 
is thus proven. 

We emphasize, however, that if the dimensionality of the system is small, the above 
condition can not be satisfied]! as 

5 > l/ds. (19) 
In particular, for a single spin 1/2 considered as a system one at best obtains from (fT^j) 

7~(7)s \ \) s <l/V2, (20) 

which is not very restrictive. Thus in the case of small d$ Theorem does not answer the 
question whether the system ISI property holds or not. Evidently in this case the system 
ISI problem requires some additional treatment. In the following section we derive a 
necessary condition for the system initial state independence, which in particular appears 
to be useful when ds is small. 



4 Necessary condition for the system initial state in- 
dependence 

First let us refine the definition of the system initial state independence. 

Definition. The equilibrium state p s of the system is independent from the initial state 

3 After the present article was completed and submitted to arXiv, we learned about a very recent 
previous work by Christian Gogolin [9] , in which he expressed the same criticism concerning the uselessness 
of the result of [B] in the case of small ds- Moreover, he proved another sufficient conditions for the system 
ISI, which works well for small ds, but relies on the eigenstate thermalization hypothesis. The latter is 
discussed in the following section. 
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of the system i/j for a fixed initial state of the bath $ with the accuracy e if 

p s [ip&] - (p s ) s [$] <e for any ip. (21) 

We remind that (...)s denotes the averaging over the normalized states from S with a 
uniform measure, while brackets [...] indicate the functional dependence. 

According to the above definition, to prove the system ISI property means to establish 
the inequality (12T1) for (almost) any initial bath state $ with some small e under reasonable 
conditions. In the present paper we do not provide such a proof. Rather we average the 
inequality (|2"TI) over $ from some subset of B and obtain a less restrictive but more 
tractable bound, which constitutes 

Theorem 1 (The necessary condition for the system ISI). Let the hamiltonian 
H have a nondegenerate energy spectrum. Let J 7 be some (possibly small) subset of 
a restricted subspace Br C B. Assume that the equilibrium state p s of the system is 
independent from the initial state of the system with the accuracy e for all initial states 
of the bath which belong to J 7 . Then 



sup 



(p s )bM-(p S )s®b r <e' (22) 



with 



e'<e + 2J^ + -±= + -e- c ^, c=1/(18tt 3 ). (23) 
V d R yd R v 

Here p is the measure of J 7 (with respect to the uniform normalized measure on Br). In 
other words, for a random $ G Br 

p = Pr {$GJ}. 

Note, however, that the dimensionality of J 7 C Br should be equal to the dimensionality 
of Br, otherwise p = and e' — oo. 

The proof of Theorem 1, which is largely based on Theorem 0, may be found in the 
Appendix. 

According to (1231) . the subset J 7 may have an exponentially small measure p, and still 
e' would be small enough to make the bound (1221) restrictive. Indeed, e' is small as long 
as £ is small and oIr is sufficiently large to ensure that p ^> e~ cy ^. Thus Theorem 1 
states that if the system ISI property holds for at least exponentially small number of the 
bath initial conditions $, then the restrictive bound (1221) is valid. Usually it is natural 
to demand that the the system ISI property holds for those initial states of environment 
which have well-defined energy. In this case the set J 7 may be constructed from those $ 
which provide a small dispersion to H B . 

The averages (p s )b r W\ an d (p s )s®B R take more explicit form in the specific case when 
B R = B. 
Lemma. 

(7)s®B = (dsy 1 ! 3 . (24) 
If further the hamiltonian H has a nondegenerate energy spectrum, then 

S*>bM = T (25) 

UK — 
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With this Lemma in hand one may reformulate Theorem 1 to obtain the following. 
Theorem V . Let the hamiltonian H have a nondegenerate energy spectrum. Assume 
that the equilibrium state p s of the system is independent from the initial state of the 
system with the accuracy e for all initial states of the bath $ from some (possibly small) 
subset T C B. Then 



sup 

V>es 



< e', (26) 



where e' is bounded according to (123]) with cIr = d&. 

Proof. Theorem 1' follows directly from Theorem 1 and Lemma. Therefore it is sufficient 
to prove the Lemma. As far as p 5 [■?/>$] is a quadratic form with respect both to ip and 
to $, the averaging over S and B with a uniform measure is equivalent to the averaging 
over arbitrary orthonormal bases in S and B correspondingly: 

ds drs 
3=1 1=1 

Applying this rule to the decomposition f[TT]) and taking into account that 

^El^nlWI^^^IPnl^) (28) 
1=1 

one obtains the equalities ( 12^|) and (125]) . □ 

Although Theorem 1 is stronger than Theorem 1', the latter may be easier applied for 
the analysis of the specific models. For this reason we concentrate on Theorem 1' in what 
follows. 

In fact Theorem 1' states that if the system ISI property holds, then the majority of 
= tr J g|^ f n )(^ r n | should be approximately proportional to the unit matrix, 

Pn * d s H s . (29) 

This requirement is natural. Indeed, according to [5] almost all vectors \l/ from "H yield 
trgl^} ($| ~ cTc, 1 ! 5 . More precisely, for a random vector \1> G "H 

Pr||| tr B |*>(*| ~d s l t s || > + e | < 2e~ cdBe2 . (30) 

Therefore, for a generic Hamiltonian H one expects the bound (126]) to hold with a fairly 
small e. 

To get more insight in the statement of Theorem V let us consider a situation in which 
the system ISI property is known to hold. Namely, consider the weak interaction case and 
assume that thermalization occurs at the level of individual eigenstates [H][l2], which 
means that pf t = Z~ l exp(— (5 n H s ) for (almost) all n% In fact general considerations and 



The inverse temperature /3„ for individual eigenstates ty n of the composite system is defined in a 

l In r I 

dE \E=E„ 



usual way, j3 n = ^f\ F=F , where the state density function r(E) is reasonably smoothed 
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numerical studies suggest that this eigenstate thermalization hypothesis holds generically, 
see e.g. [IB] . In this case all four thermalization properties are valid. In particular, ac- 
cording to eq.f llip the equilibrium state of the system is of the Boltzmann-Gibbs canonical 
form and does not depend on the initial state of the system, provided the initial state of 
the composite system has a small energy dispersion^ Let us make sure that our necessary 
condition of the system ISI holds in this case. Note that as far as the dimensionality of 
% is finite, negative temperatures are allowed as well as positive (see [H] for the discus- 
sion of statistical physics with negative temperatures). Normally in such situation the 
inverse temperature is close to zero for the vast majority of states. This is especially the 
case when the bath is composed of many weakly interacting subsystems with identical 
spectrum, as may be shown with the use of the central limit theorem. Thus the major 
contribution to the average over B comes from the states with high temperature, (3 n ~ 0. 
As a results, eq. (1291) is satisfied for the majority of n and the statement of Theorem 1' 
holds with some small e. 

Here we would like to make the following remark. Although we assume that the 
dimensionality of the Hilbert space of the bath is finite throughout the present paper, it 
seems plausible that all our results, in particular, Theorems 1, 1' (and also Theorem 2, see 
below) may be generalized to the case when alg = oo. Indeed, does not enter Theorem 

d s d B 

1 at all, while it enters Theorem 1' only through the average (c^g) -1 ^2 (iplPnl^Pni which 

n=l 

presumably remains well-defined when oIq — > oo. In this case Theorem 1' in fact provides a 
necessary condition for the system ISI in the hight-temperature regime, when it is natural 
to expect that all eigenstates of p s are equiprobable independently of the initial state of 
the system. 

Although our necessary condition for the system ISI is expected to hold for generic 
Hamiltonians H, as is clear from the above discussion, it does not hold for some specific 
Hamiltonians. This is exemplified in the next section. A simple model is discussed there 
when this condition turns out to be restrictive enough to prove the absence of the system 
ISI, although the decoherence and the equilibration occur and the bath ISI is present. 

Before we turn to the specific example let us reformulate our general results in the 
extreme case when the system S is represented by a single spin 1/2. In this case the 
equality (125]) may be rewritten as 



n=l 



(P)S = d 1 ^Pn(P0,Pn), (31) 

while the inequality f[26l) - as 

G^^PntPCbPn) 



sup |(p)e| = sup 

Po PC) 



n=l 



< ef. (32) 



It turns out that one may get rid of supremum in (I3"2"|) and obtain the following. 



5 In ref. [5] (Theorem 2.8.1) one may find a quantitative bound on the time-averaged distance between 
two states of a system corresponding to two different initial states. This bound is restrictive whenever 
the eigenstate thermalization hypothesis is valid. 
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Theorem 2. Consider the hamiltonian H with the nondegenerate energy spectrum. 
Assume that the equilibrium state of the spin S is independent from the initial state of 
the spin p with the accuracy e for all initial states of the bath $ from some (possibly 
small) subset T C B. Then 

(i) 

' 1 d d \ V 2 

^^^(P«,p m ) 2 <V3e' (33) 

v n=l m=l / 

(") 

^Ep« <3£ ' ( 34 ) 

n=l 

with e' bounded according to fl23|) with = fig. 

The proof of Theorem 2 may be found in the Appendix. 

The second bound in Theorem 2 is weaker but more tractable than the first one. It 
shows that the purities of eigenstates, tr(p^) 2 = (1 + p^)/2, should be on average very 
close to its minimal value 1/2. This requirement also enters the sufficient condition for 
the system ISI, cf. eq.fH 



5 Specific model 

In this section we concentrate on a specific class of exactly solvable (to some extent) 
models in which the system is represented by the spin 1/2 and the above derived necessary 
condition for the system ISI is not fulfilled. We consider the Hamiltonian 

H =\°* + \ J2 °«V a + H B , (35) 

a=x,y,z 

where o~ Q acts in S, V a acts in £>, at least one of V XjV is nontrivial (i.e. not zero and not 
proportional to the unit operator), all V a commute with each other, 

[V a ,V fi } = 0, (36) 

and the interaction Hamiltonian commutes with the bath self-Hamiltonian, 

[14, H B ] = Va. (37) 

Note, however, that the interaction Hamiltonian does not commute with the system self- 
Hamiltonian: [H S 1 H SB ] = (i/2)(a y V x — o~ x V y ) ^ 0. This means in particular that the 
system energy is not a conserved quantity, and one may expect some sort of thermalization 
of the system. 

Let I — 1, ds be the common eigenvectors of V a , a — x,y, z, and H B : 

V a ^ = v la $ h H B ^ = Ef^. (38) 
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The eigenvectors and eigenvalues of the total Hamiltonian read 



tfj± = V*±*i, E l± = Ef ±2^/(w + vi z y + vf x + vf y , (39) 

where ipi± are two eigenvectors of the I— dependent matrix (w<r 2 + via). We assume that 
the total Hamiltonian has nondegenerate energy gaps (and, consequently, nondegenerate 
spectrum), which is clearly a generic case. 

A specific version of the model under consideration (with V y = V z = and the bath 
composed of noninteracting spins) was introduced in [TB] in the context of decoherence 
studies. It was shown in [16] that the decoherence occurs effectively in the sense that 
the spin which is initially in a pure state rapidly becomes entangled with the bath in the 
course of the evolution. 

Equilibration and the bath ISI are also present in the model for almost all initial states 
of the bath, which follows from the general results of [B] . Namely, let us fix the initial state 
of the system ip and choose some initial state of the bath $ from a large bath subspace 
Br C B. As was proven in [6J the time averaged distance between p s (t) and p s is small, 

W^) z 7\\ < 2-^= = 4= ( 4 °) 

for almost all $ G Br. The exceptional which violate the above bound, form a set 
of exponentially small measure. This proves the equilibration property. The bath ISI 
property may be proven with the use of Theorem 0, see section [21 

Thus our model is nontrivial enough to decohere and equilibrate effectively and to have 
equilibrium states which are almost independent from the bath initial states. However, it 
does not match the necessary condition for the system initial state independence imposed 
by Theorem 2. This is essentially because the eigenstates of the composite system are 
factorized, see eq. fl3"9~|) . which results in |p n | = 1 for any n. This ensures that if e' < | 
then the restriction (l34"j) imposed by Theorem 2 can not be satisfied. In other words, for 
the overwhelming majority of the initial states of the bath the equilibrium state p s of the 
system can not be independent from the initial state of the system ip with the accuracy 
better than 1/3. 

Two other examples when the system ISI property is absent were already considered in 
[6j. In the first considered case there exist at least one conserved quantity of the system, 
i.e. a nontrivial operator A which acts in S and commutes with the total Hamiltonian H. 
Evidently equilibrium states are different for different expectation values [ip\A\ip). 

In the second case the range of energies of the self-Hamiltonian H s is greater than the 
range of energies of the combined interaction-bath Hamiltonian H SB + H B : 

F s _ F s > t?SB+B _ pSB+B 
max min "max min 

In this case the system can not transfer to (or from) the bath a substantial amount of 
energy, and the equilibrium state depends on the initial energy of the system (although 
the energy of the system is not strictly conserved). 

Our example differs from the examples provided in [6 J . Indeed, there are no conserved 
quantity of the system in our model, and the range of energies of the self-Hamiltonian H s 



12 



(which is equal to ui) may be arbitrary small. The distinctive feature of the considered 
model, which leads to the absence of the system ISI, is total lack of the entanglement of 
eigenvectors of H. Remind that the high degree of entanglement was required to prove 
the system ISI property in case when ds ^> 1 [6] (see the discussion in section [3]). Now 
we show that the absence of entanglement leads to the breakdown of the system initial 
state independence in the opposite case when dg = 2. Thus the entanglement seems to 
be an indispensable condition for the system ISI. 

We emphasize however that the exact lack of entanglement (|p n | = 1) is not of key im- 
portance in the above considerations which proved the absence of the system ISI property. 
Rather, according to Theorem 2 (ii), the value of the average ^ Ylt=i Pn * s essential. If it 
is greater than some x, then the equilibrium state of the system can not be independent 
from the initial state with the accuracy considerably better than x/3. 

6 Summary 

To conclude, we have considered the system initial state independence property - one of 
the cornerstones of the equilibrium statistical mechanics. We present a necessary condition 
for this property to hold (Theorem 1). This condition may be applied in particular in 
the case when ds (the dimensionality of the Hilbert space of the system which undergoes 
thermal relaxation) is small. This case is of special interest, as the sufficient condition 
proved previously [6] does not work for small ds- 

If we demand that the system ISI property holds with a fixed accuracy for the whole 
range of the bath "macrostates" (e.g., for all states of the bath with small energy dis- 
persion), then we get a more explicit form of the necessary condition (Theorem 1'). The 
latter indicates that the majority of eigenstates of the total Hamiltonian (which includes 
self-Hamiltonians of the system and the bath, as well as the interaction term) should be 
highly entangled. 

When the equilibrating system is just a single spin 1/2, our necessary condition leads 
to the transparent bounds on the polarization vectors of the total Hamiltonian eigenstates 
(Theorem 2). The usefulness of the derived bounds is demonstrated in the specific case. 
Namely, it is shown that for a specific form of interaction between the spin and the bath 
the necessary condition is not satisfied and thus the system initial state independence 
property does not hold. The considered interaction is not completely trivial; in particular 
it leads to the decoherence of the spin. Two other properties which are associated with 
thermal relaxation - the equilibration and the bath initial state independence - also hold 
in the considered model. 

Our results are negative in the sense that they allow only to pinpoint those models 
which lack the system initial state independence property. Further work is necessary to 
obtain more insight in the problem, in particular, to find an efficient sufficient condition 
for the system initial state independence in the case when ds is small. Also it is desirable 
to accurately generalize the obtained results to the case when the Hilbert space of the 
bath is infinitely-dimensional. 
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Appendix 

Proof of Theorem 1. We need to derive the bound (1221) from the following inequality: 



<e \fijj G S, V$ G J 7 . 



The latter along with the triangle inequality implies that 



(p S ) 



((p s )s) 



(41) 



(42) 



Now we have to move from averaging over small subset J 7 C Br to averaging over the 
whole large Br. From (T42|) one gets 



(P s )i 



((P s )s)i 



< £- 



(P S ) 



(P s )em 



+ 



((p s )s)t-((p s )s)i 



Now two last terms in the r.h.s. should be bounded. First we note that 



((P s ) S )r-((P s )s)i 



< 



< 



p s m _ {pS)l 

7m - (7h 



Tl S 



W e S. 

(43) 

(44) 
(45) 



Next we fix some tp, take some arbitrary e > and divide the set T in two nonintersecting 
parts, T\ and Ti, such as 



T 2 = < * G J 7 : 
According to Theorem (ii 



m(JF 2 ) < 2e 



> 



-cd R e 2 




(46) 



(47) 



where m(A) is the measure of the set A (remind that we take m(B^) = 1). Evidently, 



(...) 



T 



m(J'i) m(.F 2 ) 

FTV- M H 71=rV--/^2- 



and one gets 

- (?)»[ 



< 



p s m _ {pS)l 



m(JF) 



.Fi 



< < + e + ±e~ cdRe2 

dR P 



(48) 

7\m - (7) bI^p] 

(49) 
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for any ip and e, where the definition p = m(J r ) is taken into account. Now we have to 
choose the optimal one. Inserting the estimate (l4"5]) into eqs. (j4"4"|) . (l4"5"|) one evaluates the 
r.h.s. of the inequality ( ]4"3"]1 and gets the desired final expression for e' . □ 

Proof of Theorem 2. 

(i) To get the bound f l33|) from eq. fl32|) one needs to prove that 

^ d d 

SUpi^Po^PO^PO*) > - ^2^2(Pn,P m ) 2 , (50) 
P° 6 n=l m=l 

where K(p 0x ,p 0y ,p 0z ) = Yl=i E m =i(P^ Po)(Pm, Po)(Pn, Pm) = ^ 2 (p)I is a positive 
semidefinite quadratic form. One may rotate the basis to make K diagonal: 

K(P0x', POy'j POz') = ^x'Pox' + "VPo?/ + <VP(L'> — ^s 7 — <V — ■ (51) 

The maximal value of K on the unit sphere is \ z i > (A a ./ + Ay + A z /)/3 = tri^/3 = 
I Z)n=i E m =i(P'!' Pm) 2 , which is exactly the bound ([50]). 

(ii) To derive the bound ( I34p from the bound (1531) one needs to prove that 

EDp-P™) 2 ^(Ep«) • ( 52 ) 

n=l m=l \n=l / 

Let us consider the l.h.s. of the above inequality as a function of 3d variables pi, ...pd 
and find its minimum subject to d constraints of the form = a n , n = 1, ...,d, where 
< a n < 1 are some fixed numbers. We introducing d Lagrange multipliers rj n to get 
the Lagrange function L(pi, p d , T] h rj d ) = E2.=i El=i(Pn, Pm) 2 + ELi^(Pn ~ 
a n ). Differentiation of the latter over p ma gives 3d equations which (along with the d 
constraints) define the critical points: 

d 

2 ^(Pn, Pm)Pna = VmPma, TYl = 1, d, CH = X, y, Z. (53) 
n=l 

Assume that we already know the set of vectors p n which minimize the l.h.s. of eq. (|52j) 
subject to the imposed constraints. This set of vectors should obey equations (153")) . which 
may be rewritten as 

2j M a/3 p m/3 = rj m p ma (54) 

fi=x,y,z 

for every m and a. Here | \M a p\ | = ||2 Y^n=i PnaPn^l I is a 3 x 3 real symmetric matrix. Note 
that it does not depend on m, which is of key importance for the present proof. It has 
three orthonormal eigenvectors £1,^2, £3- According to (154"]) every nonzero p n is collinear 
to one of this eigenvectors and, consequently, orthogonal to two other eigenvectors. In 
other words, in the set of vectors p n which minimize the l.h.s. of eq. (l52p subject to the 
constraints each two vectors are either collinear, or orthogonal. Without loss of generality 
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we assume that pi,...,p dl are collinear with £x, p dl+1 , p d2 are collinear with and 
Pd 2 +i, •••) Pd ar e collinear with £ 3 . Then 

d d / dx \ 2 / d 2 \ 2 / d \ 2 I f d \ 2 

^^(Pn,Pm) 2 = I J^fln J +( ) + ( a ™ ) - 3 ( G " ) ' ( 55 ) 

n=l m=l \n=l / \n=<fi+l / \n=d 2 + l J \n=l J 

which proves the inequality (1521 . □ 
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